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Motivated by the potential chiral spin liquid in the metallic spin ice Pr2lr2 07, we consider how 
such a chiral state might be selected from the spin ice manifold. We propose that chiral fluctuations 
of the conducting Ir moments promote ferro-chiral couplings between the local Pr moments, as a 
chiral analogue of the magnetic RKKY effect. Pr2lr2 07 provides an ideal setting to explore such 
a chiral RKKY effect^ given the inherent chirality of the spin-ice manifold. We use a slave-rotor 
calculation on the pyrochlore lattice to estimate the sign and magnitude of the chiral coupling, and 
find it can easily explain the 1.5K transition to a ferro-chiral state. 



I. INTRODUCTION 

Pr2lr2 07 sits at the intersection of two recent fields of 
high interest in condensed matter physics: spin ice and 
iridate physic^^. Pr and Ir sit on interpenetrating py- 
rochlore lattices of corner sharing tetrahedra, where the 
Pr ions are Ising local moments while the Ir form a corre- 
lated conduction band. In such a heavy fermion material, 
one expects magnetic order or a heavy Fermi liquid, but 
there is no sign of either. Instead, below 1.5K there is an 
anomalous Hall effect that persists in the absence of both 
magnetic field and any observable magnetization. This 
hysteretic signal implies the onset of an unknown phase, 
originally proposed as a chiral spin liquid^. 

Pr^+ is a 4f^ ion with a J = 4 magnetic non-Kramers 
doublet iPsi) ^\Jz = ±4) ground stat^. The Pr mo- 
ments are therefore perfectly Ising and point into and out 
of the tetrahedra along the local (111) axis, while the in- 
plane components of the doublet are strictly quadrupo- 
lar. The RKKY interaction gives a ferromagnetic Pr- 
Pr coupling, which means the Pr sublattice realizes the 
spin ice Hamiltonian. Spin ice is a degenerate ground 
state on the pyrochlore lattices containing all possible 
arrangements of "2 in - 2 out" tetrahedra^ (see Fig. lA). 
Indeed, spin-ice correlations are seen in the magnetiza- 
tion, where a meta-magnetic transition to the "3 in - 1 
out" state is seen for fields along the [111] direction, but 
not [100] and [110], where the field simply aligns "2 in - 
2 out" tetrahedra^. The resistivity minimum at 25K is 
also consistent with the development of spin-ice^. The 
data is strongly suggestive of spin-ice correlations, but 
not conclusive in the way pinch-points in neutron scat- 
tering would bd^. However, spin-ice is clearly not the 
whole story, as Pr2lr207 has no extensive ground state 
entropy, implying further correlations. 

The Ir^+ configuration is 5d^, equivalent to one hole 
in the i2g orbital. Strong spin-orbit coupling is predicted 
to convert this configuration to a half- filled Jeff = 1/2 
doublet^. As there are 4 Ir/unit cell, semi- metallic 
behavior is expected. However, due to the narrowed 
Jeff = 1/2 band- width, the smaller 5d U is capable of 
inducing Mott insulating behavior, as seen in Sr2lr04^ 
and Na2lr03^. Pr2lr207 is part of a family of pyrochlore 
iridates, R2lr2 07 that undergo a metal-insulator transi- 
tion (MIT) as a function of the rare-earth elemenlPH^. 



Numerical studies of the pyrochlore iridates predict not 
only metallic and Mott insulating states, but an inter- 
mediary Weyl semi-metal^^. While R=Pr is metallic, 
R=Nd-Ho have MITs at temperatures increasing from 
36K to 150K. Pressure experiments on Nd2lr207^^ and 
Eu2lr207'^ have been able to suppress the insulating 
phase, indicating that Pr2lr2 07 is quite close to a MIT 
induced by increasing interaction strengths. The exact 
nature of this transition is unclear; the second order na- 
ture suggests that it has a large Slater component, but 
analogy with other iridates suggests that it likely has 
some Mott character. 

The Pr moments and Ir electrons are coupled by the 
Kondo effect, but the non-Kramers nature means that 
valence fluctuations are to excited Kramers doublets^ gen- 
erating a two-channel Kondo effect instead of the usual 
one-channel effect. Here, the symmetry of the two chan- 
nels is protected by time-reversal symmetry. While the 
ordinary Kondo effect leads to heavy Fermi liquid forma- 
tion, the two-channel version is critical and c annot form a 
Fermi liquid without a phase transition^'I^. The diverg- 
ing Sommerfeld coefficient, 7(T) = C(T)/T, the partial 
quenching of the Pr moment seen in the susceptibility 
and the resistivity minimum are all consistent with two- 
channel Kondo physics^^ While two-channel Kondo 
physics can lead to broken symmetry states, as proposed 
for URu2Si2[17 , the main signature of these hastatic 
states is the development of a heavy Fermi liquid at a 
mean- field-like phase transition. In Pr2lr207, x 7 
continue to diverge and there is no sharp phase transi- 
tion in the thermodynamic quantities. 

The absence of thermodynamic anomalies is also the 
strongest evidence against magnetic order. It is sup- 
ported by the absence of spin-freezing in /iSR datsP^, 
although, as the muons modify the Pr environment, neu- 
tron or x-ray experiments would be more conclusive. 

There is, however, a clear signature in the Hall ef- 
fect, which shows hysteresis in field and develops a 
zero- field contribution sharply at 1.5K, while the sus- 
ceptibility does not develop hysteresis until much lower 
temperatureP*. The experimentalists proposed that this 
large Hall effect is due to a macroscopic scalar spin chiral- 
ity. The scalar spin chirality, Kijk = Si-SjxSk is the solid 
angle between three spins on a triangular plaquette - as 
such it is only nonzero when the three spins are neither 
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collinear nor coplanar. In metallic systems, the chirality's 
identity as a solid angle means that conduction electrons 
coupled to the local moments acquire a Berry phase as 
they move through the lattice. This Berry phase acts as 
a local magnetic flux and contributes to the anomalous 
Hall effect, pxy = RqB -\- RmM -\- R^k^ even in the ab- 
sence of both magnetic field and magnetization^^. This 
chiral contribution can develop in magnetically ordered 
states, like Nd2Mo207[42 , but does not require magnetic 
order. Several recent theories for Pr2lr207 have pro posed 
that Rm is anomalously large, while k, is irrelevantP^^^. 
However, as 5 = and m is unobservably small, we take 
the experimentalists' point of view that the chirality is 
the main event. This picture can be checked - the Kerr 
effect can in principle measure the broken time-reversal, 
and resonant inelastic x-ray scattering (RIXS) should be 
able to detect spin-chirality directljP^. 

The Ising pyrochlore lattice has a lot of inherent spin 
chirality due to the non-coplanar Ising axes, which give 
every triangular plaquette a non-zero chirality. But is 
it possible to realize chirality without magnetization? 
While the chirality of the tetrahedron is always propor- 
tional to its magnetization, we can construct superposi- 
tions of states with large chirality and zero magnetization 
where the chirality originates from the hexagons of the 
kagome planes. However, a simple symmetry analysis 
actually reveals that any chiral state will automatically 
also be magnetic. 

In three dimensions, the scalar chirality acquires a di- 
rection given by the normal to the plaquette. The point 
group of the pyrochlore lattice is Oh, and both chirality 
and magnetization transform as the three-dimensional 
r4 irrep. Therefore a linear coupling between the two, 
7m • is expected. This problem is peculiar to the py- 
rochlore lattice, as generally the chirality and magneti- 
zation transform differently. One might expect that a 
different time-reversal symmetry breaking order param- 
eter, for example, one with symmetry could explain 
the lack of magnetization, however the existence of an 
anomalous Hall effect requires a r4 order parameter, as 
the conduction electron Berry phase transforms as 
Experimentally, it appears that there is a large chiral- 
ity without a large magnetization, so it is an important 
open question why the symmetry allowed linear coupling 
appears to be so small in Pr2lr207. 

In a 2D system, the chirality is an Ising quantity, and 
as such a chirality transition should be in the Ising (0(1)) 
universality class, with a diverging specific heat. How- 
ever, in cubic symmetry, the chirality can point along 
eight symmetry equivalent directions, meaning that it 
should instead be described as an 0(3) vector order pa- 
rameter with strong pinning to those directions, which 
gives a kink in the specific heat. The specific heat is 
presumably broadened by disorder in the real material, 
leading to the observed broad peak. 

As the spin-ice manifold contains chiral as well as non- 
chiral states, we propose that some interaction favors 
these chiral states, leading to the development of a cor- 




FIG. 1. (Color online) (A) For any given tetrahedra, the net 
chirality (a vector sum over the faces) will be proportional to 
the net magnetization. (B) Thus for the system to exhibit net 
magnetization without net chirality, the chirality must orig- 
inate from independent triangles, like the isosceles triangles 
that make up the hexagons of the kagome planes, whose chi- 
rality will have an associated direction given by the normal 
to the triangular plaquette. 



related state with a large chirality. This state is essen- 
tially a chiral spin liquid, except for the small parasitic 
magnetic moment guaranteed by symmetry - we will call 
this state a CSL*. A chiral spin liquid is a magnetic 
state that breaks time-reversal symmetry, but no oth- 
ers; specifically, it has no long range magnetic oide^^^ 
but does have a uniform (hz). Theorists originally dis- 
cussed quantum chiral spin liquids^^ , but it is also 
possible to have finite temperature classical chiral spin 
liquids, which would occur above low temperature mag- 
netically ordered chiral states. The initial theoretical 
proposals were in insulating two-dimensional Heisenberg 
magnetJ^, so it is quite interesting to find a strongly re- 
lated state in a metallic three-dimensional Ising magnet. 
This CSL* state is not differentiated from a ferromag- 
netic state by symmetry, although a topological phase 
transition is possible. In order to understand the CSL*, 
we must first understand its Hamiltonian. The natural 
question to ask is: what kind of perturbation can one add 
to the magnetic Hamiltonian to select this chiral state 
out of the spin-ice manifold? There are two generic pos- 
sibilities for Pr2lr2 07: the RKKY interaction will gen- 
erate further neighbor interactions or quantum exchange 
terms, which are enhanced over those in canonical spin 
ice by both the smaller Pr moment and the Kondo effect, 
can kinetically favor certain configurations. 

The magnetic interactions between Pr ions originates 
from the RKKY interaction^^, where the non-Kramers 
nature of the Pr moments ensures that the magnetic 
RKKY interaction is purely Ising-like, as the in-plane 



moment of the Ts doublet is quadrupolar in nature; we 
are neglecting the quadrupolar exchange. The RKKY 
interaction is long-ranged, decaying as just like the 
dipolar interaction. In insulating spin ice, the dipolar in- 
teractions delicately balance to preserve the degenerate 
spin ice state"^, however, the RKKY interaction does not 
generically do s(P. While the further neighbor interac- 
tions can select states out of the spin-ice manifold, there 
is no natural way to select a non-magnetic chiral state 
- even magnetically ordered states with chirality, but no 
magnetization require large unit cells. Furthermore, in 
three dimensions we do not expect the magnetic and chi- 
ral ordering temperatures to split. So further neighbor 
RKKY interactions do not provide a satisfying explana- 
tion for a CSL*. 

What about non-Ising exchange terms, which have 
been predicted to generate spin liquid states^^? The 
above picture is purely classical, which works well for 
the large moments in insulating spin-ice, /ieff ^ 10/iB, 
but here the moments, /ieff ^ ^fJ^B are much smaller, and 
quantum effects likely become important. Quantum ef- 
fects can be treated by perturbatively introducing quan- 
tum exchange terms, Xl(ij) J±Si-^Sj- + H.c; the specific 
origin is addressed by Chen and Hermel^^. This prob- 
lem was originally considered by Hermele, Fisher and 
Balent^, who suggested such a Hamiltonian might lead 
to a U{1) quantum spin liquid. The main idea is that 
degenerate perturbation theory favors states that can be 
connected by fiipping every spin in a loop. The first such 
term comes from third order perturbation theory and 
involves fiipping every spin in a hexagon, favoring "fiip- 
pable hexagons," where "in" and "out" spins alternate, 
as shown in Fig 2. As these hexagons have zero net chi- 
rality, the U{1) spin liquid is not a chiral spin liqui d, al- 
though it may be related to materials like Yb2Ti2 vl^^ * ^"'^l 
By inspection, the sixth order loops also do not favor chi- 
ral configurations, so the chiral ordering temperature is 
at best 0{J\/ {2Jz)^). Even taking = Jz^ this term is 
.0047^, three orders of magnitude too small. So quantum 
fiuctuations cannot generate a chiral spin liquid state at 
Th = 1.5K. 

In the next section, we explain how local moment chi- 
rality interactions can be mediated by chirality fiuctua- 
tions of the conduction electrons, giving rise to a chiral 
RKKY effect in section II. Then in section III, we esti- 
mate the magnitude and sign of the chiral RKKY effect 
for Ir electrons on the pyrochlore lattice. Finally, in sec- 
tion IV, we discuss the properties of this kind of chiral 
spin liquid and how it may be tested. 



II. CHIRAL RKKY EFFECT 

As magnetic interactions are unlikely to generate the 
Th = 1.5K transition in Pr2lr207, we must look else- 
where for the origin of the chiral state: to the conduc- 
tion electrons. Just as conduction electrons carrying 
spin mediate the magnetic RKKY effect between local 
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FIG. 2. (Color online) Hexagons which contain alternating 
"in" and "out" spins contain local modes in which every spin 
around the hexagon is flipped, resulting in another state in the 
spin-ice manifold. These "flippable" hexagons have zero chi- 
rality. Here, we indicate "in" spins with blue (darker) sphere 
and "out" with white (lighter) spheres, and have defined "in" 
and "out" self-consistently around the hexagon by picking the 
axes pointing into the alternating tetrahedra pointing out of 
the plane. 

moments^^, interactions between the local moment chi- 
ralities can be mediated by chirality fiuctuations of the 
conduction electrons, /^c(r). Generically, this effect will 
be captured by the Hamiltonian, 

where /^c(r) represents the conduction electron chirality, 
tzf{j) represents the local moment chirality of a plaquette 
centered at j, Xk{'^ ~ conduction electron chiral 

susceptibility, and Jcf is the coupling between the local 
moment and conduction electron chirality that originates 
from the Kondo interaction, Jcf ex J^/e^. The direct 
local moment chiral coupling, ^^jJffii ~ j)f^f{'^)f^f{j) 
is likely to be too small, however we shall show that a 
large coupling can be generated by integrating out the 
conduction electrons: 

J//(r-r') = -|Jc/|V(r-r'), (2) 

in exact analogy with the magnetic RKKY effect. Hence, 
we dub it the chiral RKKY effect. Such a term will 
be generated directly from integrating out the conduc- 
tion electrons directly from the Kondo Hamiltonian, as 
a sixth order term of order J%/Cp] however, we pre- 
fer this second order picture, both for the better un- 
derstanding of the physical mechanism and also for the 
relative simplicity of the calculation. The magnitude of 
this chiral RKKY effect can be roughly estimated from 
experimental parameters. If we take the single band 
estimate of n = 4.12 x 10^^cm~^[2] and the bare elec- 
tron mass, the Fermi energy is ep ^ lO^K. Then taking 
J RKKY = Ocw = 20i^r[2 , these parameters lead to a 
Kondo coupling, Jk — 31 OK and a chiral RKKY effect 
on the order of Jfc/cp ~ O.IK. This rough estimate gives 
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a coupling only an order of magnitude smaller than the 
actual ordering temperature. However, we can make a 
better estimate using equation (2), and we find a fer- 
rochiral coupling one to two orders of magnitude larger 
than this estimate. There are two straightforward ways 
to estimate this coupling, and we begin with a simple 
classical picture in the limit of infinite Kondo coupling 
to illustrate the general principle, before going on to a 
more microscopically motivated slave-rotor approach. 



III. ESTIMATE # 1: CLASSICAL MOMENTS 

A simple upper bound to the chiral coupling can be 
found by taking the limit of infinite Kondo coupling, 
Jk oo and assuming the Pr moments to be classical. 
The conduction electron spins then track the local mo- 
ments perfectly as they move around the lattice, seeing 
the local moments as a background gauge field, ^(r)^^. 
The local moment chirality is given by x w4(r), where 
^ ensures the chirality is dimensionless, and A is the area 
of a plaquette. Here, we take the conduction electrons to 
have a parabolic dispersion, Ck = k'^/2m — /i, and their 
action in the background gauge field is. 



S[c,A] = kBTY^ f ci 



2m 



Ck. (3) 



After integrating out the conduction electrons, the action 
becomes 



S[A] = [ 



^(q)nc(q)^(-q), 



(4) 



where nc(q) = l/(^(q)^(-q))^ = -(Jc(q) Jc(-q))^ 
is the transverse current-current correlator for the con- 
duction electrons, and we have already taken the static 
limit for simplicity. This action can be rewritten in terms 
of the chirality, A>:(q), 



»1 = ^ 



nc(q) . . 



(5) 



For a parabolic dispersion, the current-current correlator 
can be calculated analytically in the zero temperature 
limit. 



q|2/(ek+q) -/(ek) 



n 
m 



1 , r 



1 /I 



2g + g3)log 



4\q 



ek — Ck+q 
1 + 9 



1 



, (6) 



where we have introduced q = q/2kF for clarity. 

The chiral coupling constant, J/^(r) can now be calcu- 
lated 



Mr) 



^ 2 f\, ^ nc(q) ^, 

dqq I a cos 6/ 7; — e 



2TT^Ah 



dq 



Hc(g') sin2/ci?(3'r 



(7) 



The result is plotted in Fi gure 3, using hp = 27r/(1.218a) 
taken from the Hah effecP^, n = 4.12 x 10^^ cm"^ and 
the lattice constant, a = 10. 4A. The finite leads to os- 
cillatory behavior, so both ferrochiral and anti-ferrochiral 
interactions are possible, depending on the distance. For 
the parameters of Pr2lr207 , the chiral coupling at r = a 
is -200K. As we have taken Jk 00, this estimate is 
obviously an upper bound on the real value. 
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FIG. 3. The chiral coupling Jtt^(r) calculated in the classical, 
Jk 00 approximation for a parabolic dispersion. The y- 
axis is in Kelvin, for a possible set of parameters for Pr2lr2 07. 
For r = a = 10. 4A, the chiral coupling is -200K. Just like the 
magnetic RKKY interaction, the chiral RKKY interaction is 
oscillatory and decays rapidly. 



IV. ESTIMATE#2: SLAVE ROTORS 

To make a more accurate estimate, we return to equa- 
tion (2), and now attempt to calculate the conduction 
electron chiral susceptibility, XH:(r — r'). As it is a twelve 
point correlation function, this calculation initially ap- 
pears quite tedious. However, by introducing the slave- 
rotor mean-field approximation^^, we can greatly sim- 
plify the calculation, and in the process model the ap- 
proach to the Mott transition. As electrons approaching 
a Mott transition have increased chiral fluctuation^^, we 
expect the coupling to be enhanced. 

Slave-rotors are typically introduced to capture the 
quantum criticality associated with the Mott transition 
in the Hubbard model by splitting the conduction elec- 
tron into a neutral spinon and a charged holon, cj^ = 
//^e*^^[33]. This decoupling possesses a [/(I) gauge sym- 
metry. 



U(l) gauge symmetry : 



ft _^ ft 

Oi Oi-Q 



(8) 



which glues the spinon and holon back together to form 
the charged electron in the metallic phase. The metal- 
lic phase is captured by the uniform condensation of this 
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rotor, (e*^^), while it is uncondensed in the Mott insulat- 
ing phase and the Mott transition is a 4D XY transition 
before the coupling to the gauge field is taken into ac- 
count. The spin chirality is particularly simple in this 
approacli^, 



V X a = —Scl • Sc2 X Sc3, 



(9) 



where V x a is the lattice curl of the gauge field around 
a triangular plaquette, and the units are fixed by ^ as 
before. The chiral susceptibility can then be calculated 
from the transverse gauge propagator, (a(q)a(— q))^, 

xfN = (<(q)^^(-q)) 

= ^(Q'^ab - g„96)(a(q)a(-q))^. (10) 

As the spinon and holon add in series, the gauge propa- 
gator is. 



(a(q)a(-q))^ 



[n?(q) 



n5(q)]" 



(11) 



where n|i^(q) are the transverse current-current corre- 
lators for the spinous and rotors, respectively. 

We now calculate these correlators within a specific 
model: electrons on the half-filled Ir pyrochlore lat- 
tice, where the strong spin-orbit coupling means the 
5d^ Ir configuration is simplified to an isotropic effec- 
tive Jeff = 1/2 doublet^^. Both oxygen mediated Ir-0- 
Ir and direct Ir-Ir hopping are expected^^, and we take 
the nearest-neighbor dispersion considered by Witczak- 
Krempa et al^^, where the oxygen mediated hopping has 
magnitude to and the direct hopping is governed by two 
parameters, to- and t^^ with the specific choice to- = — to, 
t-jT = |to- This dispersion has no Fermi surface at half- 
filling, instead having a peculiar quadratic band touching 



with zero density of states at Ep. We add a small next- 
nearest-neighbor oxygen mediated hopping, t2 to convert 
the quadratic band touching into a semi-metallic disper- 
sion with small electron and hole pockets. The strong 
spin orbit coupling means that the Ir-Ir hopping is spin- 
dependent, leading to a complex hopping matrix r^*^,-^*^-^^^ , 
where a, 6 G {1 — 4} label the four sites per unit cell. The 
Hubbard model for this system is then. 



H 



u 



E(E4 



■E 

iajbaa' 



(12) 



As the nearest neighbor hopping has several components, 
hidden within r^^^}''^^\ ti = t/^ij-^ is actually just to- Here 
the conduction electron density is kept fixed to half-filling 
by /i, and the proximity to the Mott transition is tuned 
by changing the interaction strength, U . We introduce 
the slave-rotor decoupling, c\^^ = fjaa^^^'''^ where from 
now on we will work in terms of Xia = e^^^'^ . The hop- 
ping term is now a quartic term, and we introduce four 
Hubbard- Stratonovich fields to decouple it: 



(1,2) 
fl,2) 



^^'^^^'^\fLMViiUiim^ (13) 



where {ij) and {{ij)) are really a shorthand restricting 
{ia) and {jh) to be nearest or next-nearest neighbors, re- 
spectively. In other words, we are assuming the Qf b^^^^^ 
take one value, for nearest neighbors, another, for 
next-nearest neighbors, and are zero otherwise. This de- 
coupling leads to quadratic actions for both the fermions 
and bosons. 



^ iaa (ij) ,ab,aa {{ij)) ,ab,aa 

Sx = fdrY, + A.. - 1) - \t,Q^^^ ^l^^^ - '^Qx E ^^aX,b. (14) 

^ ia (ij),a6 {{ij)) .ab 



where \ia is a Lagrange multiplier enforcing the phase 
nature of the slave rotors that we will take to be uniform, 
= A, and we have rescaled U to preserve the 

atomic limit^^. Note that while the fermions still have the 
complex hopping matrix of the original conduction elec- 
trons, the bosons behave like spinless fermions. The near- 
est neighbor hopping element is then given by the trace of 
to^^^^i^"^^^ = to -\- ta -\- 1^^ = As the slave rotor tech- 
nique is fairly standard, and has already been applied to 
the pyrochlore lattice, with different dispersions^^ we 
have skipped several intermediate steps. The fermionic 



and bosonic bands are shown in Figure 4 A and B, re- 
spectively. The Green's functions are 

Gkmi^UJn) = . _ , (15) 

where Ckm are the eight fermion bands and ("/en are the 
four bosonic bands. We have allowed the bosons to con- 
dense in the lowest band with the uniform amplitude 
Xq = (Xia)^ which implies that the Lagrange multiplier 
A = 4ti(5^^ + l2t2Qx^ is fixed to the bottom of the 
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bosonic bands. We will concentrate on the mean- field 
theory of the metallic phase and quantum critical point 
at zero temperature, where Xq G [0, 1] quantifies the dis- 
tance from the quantum critical point. 




FIG. 4. (Color online) Bands on the pyrochlore lattice. (A) 
The bosonic bands plotted through the high symmetry points 
of the pyrochlore Brillouin zone, for t2 = -.03t, Q^^ = .14, 
Q^-p = .20, where t is the nearest neighbor oxygen- mediated 
hopping magnitude. (B) The fermionic bands with ta = — t, 
= 2/3t, and ts = -.03t. These are plotted for [/ = 0, 
and where /i = — .31t has been included so that Ep is at zero 
when the Ir sites are half-filled. (C) The fermionic density of 
states. (Inset) The pyrochlore Brillouin zone. 



The remaining parameters {Xq^Q^ 



(1,2) 



(5x'^\/i) can 



be determined from the six mean field equations, 



1 

tl,2 
Di,2 



U 



•'^ s , 

kn 



4(A + Ckn) 



Ed- 



(1,2) 



u 



3t 



1,2 



' kn V 



2Di,2K 

km 



~ ^km) 



(16) 



where Di = 6ti and D2 = 12^2 are determined by the 
relevant number of neighbors, and Afg is the total num- 
ber of sites. For ^2 ^ 0, these must be solved self- 
consistently. We choose ^2 = .05ti, and find Uc = 1.35ti. 



1 — a/^, while Qx'^^ = -14, .20 are independent of 



U. The U dependence of the renormalized bandwidths, 
^^''"^ - 1 - . /tt^, where /7,(i,2) 



1.5^1, 1.6ti; it is 



important to note that while the quasiparticle weight, 
Z = Xq goes to zero at the quantum critical point 
(QCP), the bandwidths decrease with increasing [/, but 
remain finite at the QCP. 

With the mean field parameters in hand, we can now 
calculate the current-current correlators, where the cubic 
symmetry implies n|i^(q) = 2np^^(q = qx). 



-01x(q)jlx(-q)) - {Tf,x), 

and the paramagnetic currents are given by. 



(17) 



jl(q) = T^^i(.L 

iuJn km 

jl(q) = T^^i(C 



V 



k+qm) fkmfk+qm 
k+qn)^kn^k+qn, (18) 



kn 



with the fermionic, = dekm/dky and bosonic, Vkn = 



dCkn/dky velocities. The diamagnetic contributions are. 



fkmf^ 



km 



iojn km 



dkl 

QJ^2 ^kn^kn- 



(19) 



The fermionic current-current correlator is shown in 
Figure 5 A for several values of Xq . lip decreases in mag- 
nitude with increasing U . 

The bosonic correlator requires a bit more care. Deep 
in the Fermi liquid, the paramagnetic term vanishes, and 
the diamagnetic term, T^^ = -(S/StiQi + 24^2(32)^0 
is the only contribution. At the QCP, we can calculate 
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FIG. 5. (Color online) (A) Fermionic transverse current- 
current correlator, Urici) on the pyrochlore lattice for vary- 
ing proximity to the Mott transition: Xq = (red, U = l-4t); 

= 0.05 (orange, U = 1.2t); Xi = .1 (yellow, U = Lit); 
Xi = 0.25 (green, U = 0.8t); Xq^ = 0.5 (blue, U = .3t); 

= .75 (purple, U = .08t); and Xq^ = 1 (gray, [7 = 0). Note 
that the distance from the Mott transition increases the ef- 
fective fermionic bandwidth. (B) Bosonic transverse current- 
current correlator, nx(q) on the pyrochlore lattice for vary- 
ing proximity to the Mott transition (same color scheme as 
above). Note that for large Xq, the correlator approaches a 
constant in ^-space, just as expected in the Fermi liquid. 



n^^(q) directly. In between, we must interpolate be- 
tween the twG^^, 



(8/3iiQi+24t2Q2)X2 



■ (20) 



The results are shown in Figure 5B for a variety of Xq, 
and generally decreases with increasing /7, like H^. 

Having the correlators, we now combine them and 
Fourier transform, 



" 8ttf 7r3 



1 



n^(q) + n3^(q) 



dqfiq) 



(21) 



BZ 



where we have defined A 



32 



as the area of the isosce- 



constant. We have also introduced q = qa and used that 
the units of IIf,x are t/a^ to introduce the dimension- 
less function /(^) = t/a^{q'^ - g^)/[n|;(q) + nj(q)]. For 
Jc/, we have simply taken J^P^^ where p is the density 
of states at the Fermi energy, and we take the estimate 
Jk from section II. We took t = lO^K roughly from the 
density functional theory calculation^. 




les triangle in Figure IB and a = 10. 4 A is the lattice 



FIG. 6. (Color online) Chiral coupling Jtt:(r) in real-space 
for several X§ (color scheme explained in figure 5). JK,{r) 
is ferro-chiral for all r and Xq, and increases in magnitude 
as the Mott transition is approached. The coupling constant 
is given in Kelvin, and even for [7 = 0, the chiral couplings 
between hexagons, r — aj \f2 and between layers, r — vTT/3a 
are both on the order of lOK. 



The chiral coupling oscillates slowly and is nearly al- 
ways negative, or ferro-chiral, due to the small carrier 
density, although for larger distances further from the 
Mott transition it can be small and antiferro-chiral due 
to the small pockets. By contrast, on the half-filled fee 
lattice, which has a large Fermi surface, the sign of the 
chiral coupling is very sensitive to the inter-atomic dis- 
tance, and other parameters. While the Mott transi- 
tion does increase the chiral coupling, it is already quite 
large for free fermions on the pyrochlore lattice, making 
the chiral RKKY coupling a plausible explanation for 
PrsIrsOr. 



V. PREDICTIONS AND SPECULATIONS 

As Pr2lr207 differs so substantially from the theoret- 
ical systems initially proposed to contain chiral spin liq- 
uids, determining the relevant microscopic Hamiltonian 
is an important problem. As magnetic interactions seem 
unlikely to generate the relatively high transition tem- 
perature, Th — 1.5K, we propose that the chiral RKKY 
effect generates a ferrochiral coupling. Using a slave rotor 
treatment of the Jeff = 1/2 pyrochlore lattice to estimate 
its magnitude and sign, we found that the interaction is 
always ferrochiral and of the necessary order of magni- 
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tude. The relevant Hamiltonian is then, 

{ij) ij 

Studying the ground and excited states of this Hamilto- 
nian is an interesting problem for future work, although 
unfortunately complicated by the sign problem. One key 
question is whether or not the observed smallness of the 
parasitic magnetic moment is consistent with the Hamil- 
tonian and how it might be tuned. 

As the chiral RKKY effect requires the presence of con- 
duction electrons, we do not expect the insulating ana- 
logues, Pr2Sn207[43 and Pr2Zr207[4T to exhibit chiral 
order, unlike other theories based on the Ts Pr doublelP^. 
Indeed, both these compounds appear to have residual 
magnetic entropy consistent with spin ice, although the 
non-chiral nature could be confirmed by RIXS^^. 

Pr2lr2 07's proximity to a metal insulator transition 
provides the opportunity to test whether or not driving 
the system closer to a Mott transition increases the chi- 
ral ordering temperature by increasing the conduction 
electron chiral susceptibility. The nature of the metal 
insulator transition is unclear; though its second order 
nature suggests it has a large Slater character, analogy 
to other iridates, like Sr2lr04^ suggests that it may have 
some Mott character as well. As pressure drives the ma- 
terials further from the MIT, it should suppress the chi- 
ral transition in Pr2lr207, while Nd doping should, in 
principle, increase it. It is also possible that Nd2lr207 
under pressure may realize a chiral state, and Hall effect 
measurements on single crystals under pressure should 
be done to check for such a chiral state. 




FIG. 7. Classical versus quantum chiral spin liquids. (A) 
A generic phase diagram. The classical CSL* (cCSL*) occurs 
above a magnetically ordered state with net chirality (CMO), 
while the quantum CSL* (qCSL*) is only revealed when quan- 
tum fluctuations kill the magnetic order. (B) A rough sketch 
of the frequency dependence of the spin structure factor to 
illustrate the difference between classical (solid line) spin liq- 
uids, whose fluctuations are mainly confined to cj < T and 
quantum (dashed line) chiral spin liquids, whose fluctuations 
have a much broader distribution at finite frequencies. 

An interesting open question is whether the observed 
CSL* is classical or quantum in nature. Classical spin liq- 
uids are incoherent thermal superpositions of states that 
occur above magnetically ordered states, while quantum 
spin liquids involve a coherent superposition of states. 
Generically, increasing quantum fluctuations can tune 
a thermal CSL* into a quantum one (see Figure 7A). 
The magnetic quantum criticality seen in the magnetic 
Gruneisen factor suggests that the ground state is not 
magnetically ordered, and thus Pr2lr207 may be a quan- 
tum CSL"*^. This question can be resolved by examin- 
ing the energy distribution of spin fluctuations, through 
the frequency dependence of the spin structure factor. A 
classical spin liquid is dominated by thermal fluctuations, 
so most fluctuations occur at <C t, while a quantum 
spin liquid will have a broader distribution. Theoreti- 
cally, quantum CSL*'s may have interesting topological 
properties: in the 2D case, chiral spin liquids are pre- 
dicted to have semionic excitations^^ and a quantized 
thermal Hall effect^^, and it is unclear how the three- 
dimensionality and small parasitic moment will affect the 
topological nature of the phase. 
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